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1 . 1  BACKGROUND 


We  proposed  to  initiate  a  program  of  research  in  real-time 
holographic  image  processing.  In  this  effort  new  device 
architectures  were  proposed,  supported  by  an  experimental 
demonstration  and  theoretical  investigation. 

The  new  applications  which  we  worked  on  were  speckle 
velocimeters,  incoherent-to-coherent  nonlinear  joint  transform 
correlators,  and  imaging  within  a  scattering  medium  using  homodyne 
detection  in  real-time  holography. 

For  developing  these  projects  in  nonlinear  holographic  image 
processing,  it  is  necessary  to  study: (a)  the  coupled  wave  theory  in 
real-time  holographic  devices (1)  ; (b)  the  scattering  effects  from 
specific  cases  of  random  moving  media(2);  and  (c)  the  effects  of 
holographic  nonlinearites  on  spatial  and  temporal  signals  using  the 
nonlinear  transform  method  (3)  .  Our  nonlinear  holographic  image 
pj^Qcessor  used  mainly  two  devices:  (a)  An  incoherent-to-coherent 
convertor  (4) ;  (b)  A  real  time  holographic  homodyne  and  heterodyne 
detector (5 , 6) . 
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1.2  FR06S21H  OBJECTIVES  2^  GOALS 

The  program's  objective  and  goals  were  to  use  real-time  holographic 
media,  such  as  photorefractives(l)  and  bacteriorhodopsin(2) ,  to 
understand  their  nonlinearities  and  then  to  use  these 
nonlinearities  to  implement  various  image  processing  algorithms.  In 
our  implementation  we  used  two  main  devices;  (a)  The  incoherent— to— 
coherent  detector;  (b)  The  real-time  holographic  homodyne  and 
heterodyne  detector. 

IKCOHERENT-TO-COHERBNT  CONVERTOR 

The  incoherent-to-coherent  converter  ( 18 )  acts  like  a  negative 
exposure  film.  For  low  intensities  it  behaves  as  a  square  law 
receiver  and  for  large  intensities  it  behaves  as  a  saturated  square 
law  receiver.  This  characteristic  may  be  used  to  implement  the 
nonl inear  j  oint  transform  correlator (7,8,9). 

HOLOGRAPHIC  HOMODYNE  AND  HETERODYNE  DETECTOR 

Mixing  any  two  rapidly-varying  fields  in  any  real-time 
holographic  medium  is  always  accompanied  with  time  integrative 
behavior.  This  time  integration  behavior  is  the  origin  of  homodyne 
and  heterodyne  detection  (10,11,12).  In  electronics  hbmodyne  and 
heterodyne  detection  is  the  basis  for  designing  many  devices  such 
as  the  time  integrative  correlator,  frequency  division 
demultiplexing,  lock-in  amplifier  and  notch  filter.  All  these 
electronic  applications  can  be  implemented  serially  as  well  as 
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spatially  using  real-time  holography.  We  will  use  the  holographic 
lock-in  amplifier (10, 12)  for  detecting  images  imbedded  in  noise  and 
we  will  use  the  homodyne  detector  in  four-wave  mixing  for  speckle 
velocimetry . 


1.3  TECHNICAL  PROGRAM  SUMMARY 
1.3  (a)  Inverse  Nonlinear  Joint  Transform  Correlator. 

Pattern  identification  using  optical  correlators  can  play  an 
important  role  in  applications,  such  as  target  recognition  and 
robotics  vision.  Two  main  approaches  have  been  used  to  design  the 
optical  correlator.  The  first  one  is  based  on  the  Vander  Lugt 
correlator  geometry,  the  second  is  based  on  the  joint  transform 
correlator  JTC  geometry. 

The  second  approach  is  more  convenient  in  the  sense  that  there 
is  no  need  to  fabricate  holographic  filters,  such  as  the  phase-only 
filter  or  the  matched  filter.  The  original  joint  transform 
correlator  required  a  square-law  receiver  in  the  Fourier  plane  (13) . 
This  was  the  basis  of  all  the  subsequent  designs  of  the  opto¬ 
electronic  JTCs(9).  In  contrast,  real-time  holographic 
implementation  of  the  JTC  was  more  like  the  vander-Lugt  in  terms  of 
its  transfer  function.  Only  recently  we  implemented  the  first 
holographic  joint  transform  correlator  based  on  a  square-law 
receiver  in  the  Fourier  plane  (7,8).  Square  law  nonlinearites  can 
be  implemented  in  various  forms.  Our  specific  implementation  will 
use  the  incoherent-to-coherent  converter (4) .  The  square-law 
nonlinearity  in  the  incoherent-to-coherent  convertor  has  a 
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saturative  behavior  like  a  negative  exposure  film.  This  fact  makes 
the  incoherent-to-coherent  convertor  an  ideal  device  for 
implementing  the  nonlinear  joint  transform  correlator.  There  are 
two  advantageous  features  of  our  correlator;  it  will  have  a  reduced 
DC  peak  and  the  joint  spectrum  saturation  can  be  reached  easily  by 
using  different  wavelengths  for  erasure  of  the  grating. 

1.3(b)  Speckle  velocimeter. 

A  speckle  velocimeter  is  a  device  of  interest  for 
metrological  applications.  In  the  past,  various  schemes  have  been 
used  for  these  purposes.  Some  of  these  schemes  used  the  algorithm 
of  zero-crossing  rates (14)  ,  others  used  measurement  of  the 
autocorrelation  ftinction;  while  others  used  lock-in  electronic 
holography.  Only  three  publications  using  real-time  holography  are 
available  in  the  literature  for  these  purposes (15, 16, 17) .  In  all 
these  techniques,  the  processing  was  in  real  time.  However,  the 
data  needs  further  computer  processing.  We  propose  an  alternative 
technique  which  is  based  on  using  the  holographic  homodyne  and 
heterodyne  detection (10, 12)  as  time  correlation  motion 
detect  ion  ( 17 ) .  In  this  technique  we  expect  to  have  a  direct 
measurement  of  the  velocity  or  the  angular  velocity  without  any 
further  processing. 

\ 

1.3(c)  Imaging  in  scattering  media 

Imaging  in  scattering  media  is  one  of  the  most  challenging 
problems  in  optical  signal  processing.  Holography  is  one  way  that 
appears  quite  promising.  In  fact,  there  is  a  variety  of  holographic 
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methods  that  may  be  used  for  imaging  in  scattering  media. 

The  earliest  approach  for  imaging  in  scattering  media  was  due  to 
Duguary  and  Mattick  (1971)  who  used  the  principle  of  the  first 
arriving  light  (18) .  In  this  technique  a  pulse  of  light  bearing  an 
image  passes  through  the  scattering  medium.  The  light  that  scatters 
least  emerges  first;  whereas,  light  that  scattered  more  emerges 
later.  Through  time  grting  it  was  possible  to  filter  the  scattered 
light  .  Both  real-time  holography  and  electronic  holography  have 
been  used  for  these  purposes.  Only  the  unscatterd  light  can 
interfere  with  the  plane  wave  and  form  a  hologram.  We  suggest  an 
alternative  technique  for  imaging  within  the  scattering  medium. 
This  technique  is  based  on  using  the  real  time  hologram  as  a 
spatial  lock-in  detector. 
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Chapter  2 


Incoherent  Erasure  Joint  Transform  Correlator 


2.1.0  Abstract 

We  propose  a  nonlinear  joint  transform  correlator  which  is  produced 
from  the  erasure  of  a  real  time  hologram  by  the  joint  transform 
spectrum.  Our  analysis  of  this  incoherent  erasure  correlator  with 
a  four  wave  mixing  read-out  arrangement  indicates  that  the 
performance  of  this  correlator  exceeds  that  of  the  phase  only 
correlator 
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2.1.1  Introduction 

We  have  recently  introduced  a  new  class  of  real-time  holographic 
nonlinear  joint  transform  correlators  (JTC)  operating  on 
soft-clipping  quadratic  transfer  functions  [1-2].  These 
correlators  tune  from  th®  matched— filter  to  the  phase— extraction 
limit  and  enhance  signal  detection  in  noisy  environments  by 
compressing  the  spectra  of  both  the  signal  and  the  noise. 

Our  first  approach  was  based  on  energy  transfer  in  a  two-beam 
coupling  JTC(TBCJTC)  [1]  using  barium  titanate.  However,  many  fast 
real-time  holographic  materials,  such  as  polymers  and  certain 
geometries  of  multiple  quantum  wells,  cannot  be  used  in  this  device 
because  they  do  not  produce  two-beam  coupling.  Our  second  approach 
was  based  on  four-wave  mixing  [2]  as  an  alternative  to  two-beam 
coupling,  utilizing  a  self-pumped  phase-con jugator  to  retroreflect 
the  joint  spectra.  This  geometry  is  not  attractive  because 
efficient  self-piomped  phase  conjugation  is  slow. 

We  therefore  propose  the  incoherent  erasure  JTC  (lEJTC)  based 
on  the  incoherent  erasure  of  real-time  holograms  [3]  and 
incoherent-to-coherent  conversion  [4].  Our  analysis  shows  that  the 
lEJTC  operates  with  a  tunable  nonlinearity  which  can  be  adjusted  to 
operate  from  the  matched  filter  to  the  phase  extraction  limit  (or 
inverse  filter) ,  in  a  fashion  similar  to  our  previous  approaches 
[1-2].  The  nonlinear  transfer  function  of  the  lEJTC  is  controlled 
by  the  joint  spectral  intensity  at  the  Fourier  plane.  The 
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classical  matched  filter  is  produced  in  the  joint  transform  low 
intensity  limit,  while  the  inverse  filter  is  produced  in  the  joint 
transform  high  intensity  limit.  Such  tunability  is  very  important 
for  detecting  signals  in  various  noise  environments.  For  example, 
the  matched  filter  has  proven  to  be  effective  for  additive  Gaussian 
noise  [5]  but  fails  to  detect  signals  in  clutter  that  can  be 
detected  by  the  phase  only  filter  [6]. 

The  main  operating  improvements  of  the  lEJTC  are  the  increased 
speed  and  reduced  noise  provided  by  the  use  of  a  second  wavelength 
for  the  erasure  (incoherent)  beam.  Using  a  shorter  wavelength  for 
the  erasure  light  enables  fast,  low  power  nonlinear  operation  into 
the  saturation  region,  and  allows  easy  separation  of  any  erasure 
light  noise  from  the  correlation  signal. 

2.1.2  Background  and  Device  Theory 

A  proposed  experimental  set-up  is  shown  in  Fig.  1:  A 
temporally  coherent  set  of  beams  of  wavelength  kc  establishes  a 
four-wave-mixing  arrangement.  Beams  1  and  4  write  a  grating  and 
beam  2  is  the  readout  which  produces  beam  3,  the  phase-conjugate 
output  beam.  The  erasure  beam  uses  temporally  incoherent  light  at 
the  input  plane  and  the  spectra  of  the  scene  S(y>^,v^)  and  the 
reference  R(v„,  v„)  )  are  focused  on  the  photorefractive  crystal. 
Here  Vj^  and  Vy  are  the  spatial  frequency  coordinates.  The 
generation  of  charge  carriers  by  the  absorption  of  the  joint 
spectrum  erases  the  existing  grating  and  this  erasure  depends  on 
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the  intensity  |R(Vx,Vy)  +  S(Vjj,Vy)j2. 

The  incoherent  erasure  beam  is  temporally  incoherent  to  the 
other  beams  and  must  have  a  at  least  partial  spatial  coherence  to 
produce  the  joint  transform.  This  beam  may  be  produced  by  a 
variety  of  optical  techniques  including  those  of  white  light 
processing.  We  consider  a  crystal  with  coupling  coefficient  y, 
thickness  L,  and  a  small  absorption  coefficient.  Before  the 
erasure  beam  is  introduced  the  phase”Con jugate  beam  amplitude  Ajq 
for  the  undepleted  pump  case  is  given  by  [ 7 ] : 


■^30  = 


yL 


Where  Aj  is  the  amplitude  of  the  jth  beam  and  is  the  sum  of  the 
intensities  of  the  coherent  beams  1,  2,  and  4.  The  incoherent  beam 
does  not  couple  with  the  existing  beams  due  to  strong  violation  of 
the  Bragg  conditions,  but  the  total  absorbed  intensity  increases. 
Considering  a  monochromatic  source  (e.g.  a  second  laser)  of 
wavelength  A,,-  ,  the  modulated  phase-conjugate  beam  amplitude  is 
given  by: 


A,(v 


x> 


{A^ADA^yl 


|■R(v^,Vy)  ^^^5(Vy,Vy)  P 
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where  I,-  is  the  incoherent  intensity  (of  the  erasure  beam)  before 
the  input  plane,  fl  is  the  transform  lens  focal  length,  and  is 
a  dimensionless  factor  which  depends  on  the  absorption  coefficients 
and  wavelengths  of  the  incoherent  and  coherent  beams  and  can 
strongly  enhance  the  incoherent  erasure. 

The  modulation  depth  which  describes  the  phase-conjugate 
beam  amplitude  modulation  by  the  spectral  energy,  is  obtained  by 
setting: 

A3+A30  [l-d(Vj^,  vp  ] 


Then  one  can  show  that: 


d{v^,Vy) 


mgE{\^,Vy) 
l+mgE{\^,  Vy) 


where  E(v„,v„)  is  the  normalized  energy  spectrum  and  me  is  the 
effective  beam  ratio.  This  equation  is  functionally  identical  to 
the  low— coupling  limit  of  the  transfer  function  of  the  four-wave 
mixing  JTC:  see  for  example  Eq.  (4)  in  reference  [2]. 

If  Eq  is  the  product  of  the  area  of  the  input  plane  and  its  DC 
component,  we  define  and  m^  as  follows: 
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F(Vjc,  Vy) 


|jgVjC,Vy)  +g(V;t,Vy)  1^ 

Eo 


m  -<7  ■*^0 


£i 

-Z^c 


In  Fig.  2  we  plot  the  modulation  depth  d(Vjj,vy)  versus  the 
product  mgg(VxfVy)  ),  which  is  the  relative  spectral  energy  incident 
on  the  crystal.  Two  distinct  patterns  are  clearly  identified:  the 
modulation  is  linear  with  small  energies,  while  for  large  energies 
it  is  independent  of  the  signal  energy.  This  soft-clipping 
saturation  commences  at  the  point  of  inflection,  where  me'E  =1. 
Since  the  crystal  responds  to  the  product  m^g,  not  just  the  energy 
E,  the  value  of  m^  substantially  affects  the  way  the 
phase-conjugate  beam  amplitude  is  modulated. 

For  a  small  beam  ratio  (<  m^  =  10'^  )  we  can  set  m^  «  1 
for  all  the  spectral  energy  values  and  the  modulation  depth  can  be 
approximated  as: 


d(Vjf,v_y)  ■ingJF(Vjf,Vy)  . 

Equation  7  agrees  with  Fig.  2  in  that  the  modulated  beam  has 
a  linear  dependence  with  small  signal  energies.  This  dependence  is 
typical  of  the  matched  filter,  whose  peaks  are  simply  the 
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correlation  between  the  reference  and  the  scene. 
For  a  large  beam  ratio  (>  me  =  10^  )  we  can  set  m^ 


and 


the  modulation 


depth  can  be  approximated  as: 


d(v  ,vj  =1- —  ^  ■■ 

*  y  jn„E(v„,vJ 


This  response  is  independent  of  the  signal  intensity  at  the  limit  me-®' 

is  equivalent  to  the  inverse  filter  and  leads  to  phase  extraction. 
As  the  beam  ratio  increases  the  nonlinearity  at  the  Fourier  plane 
is  transformed  from  the  matched-filter  type  (the  linear  section)  to 
the  inverse-filter  type  (the  limiting  section) .  For  intermediate 
beam  ratios  the  portions  of  the  joint  spectrum  with  low  energy, 
the  high  frequencies,  are  amplified  by  the  linear  region,  and  the 
ones  with  high  energy  (low  frequencies)  are  deamplified  by  the 
limiting  region. 

2.1.3  Computer  Simulation 

We  ran  computer  simulations  of  the  lEJTC  with  Ajq  =1  and  m^  set 
to  10’^,  10^,  and  10^.  We  correlated  two  disks  (Fig.  3a),  and  a 

tank  versus  a  tank  in  clutter  (Fig.  4a).  Figs.  3b-c  plot  the  peak 
intensity  arrays  produced  by  the  matched  and  phase-only  filters 
with  the  disks.  Figs.  3d-f  display  the  peaks  produced  by  the  lEJTC. 
In  Fig.  4b-f  the  corresponding  peaks  with  the  tanks  are  plotted. 

The  shape  of  the  peaks  and  the  intensity  strongly  depend  on 
the  beam  ratio,  which  verifies  the  theoretical  predictions.  The 
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correlator  with  small  m*  values  has  characteristics  of  the  matched 
filter.  We  note  the  similar  shape  of  the  correlation  peaks  between 
the  lEJTC  with  me=  10*^  and  the  matched  filter.  Both  fail  to 
produce  a  correlation  peak  with  the  tank  in  clutter. 

transformation  of  the  nonlinearity  from  the  matched~filter  type 
to  the  inverse-filter  type  is  indicated  by  peak  sharpness  and 
intensity  that  progressively  exceed  those  of  the  phase-only  filter 
as  the  beam  ratio  increases.  For  example,  the  correlation  peak  for 
the  disks  produced  by  the  lEJTC  with  me  =  10^  had  PNR  [9]  89.66 
and  intensity  308.  The  peak  produced  by  the  matched  filter  had  PNR 
4.75  and  intensity  11.3,  and  the  peak  produced  by  the  phase-only 
filter  had  PNR  22.32  and  intensity  13.5. 

2.1.4  Conclusion 

We  have  presented  a  new  approach  to  implement  the  nonlinear  joint 
transform  correlator.  The  lEJTC  is  based  on  the  erasure  of  a 
four— wave  mixing  grating  in  a  real-time  holographic  media.  The 
resulting  soft-clipping  quadratic  transfer  function  is  similar  to 
our  previous  photorefractive  implementations.  The  operation  of 

these  nonlinear  correlators  can  be  tuned  from  classical  matched  to 

/ 

the  phase  extraction  by  appropriate  selection  of  the  operating 
point.  In  addition,  low-power  operation  may  be  achieved  by  the 
proper  choice  of  the  wavelength  of  the  erasure  beam. 
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2.1.5  Figure  captions: 

(1)  Proposed  experimental  arrangement. 

(2)  Modulation  depth  of  the  phase-conjugate  beam  amplitude  versus 
the  relative  spectral  energy. 

(3)  (a)  The  input  plane  with  disks  used  for  the  simulation 

Correlation  peaks  intensity  plots  for  the  two  disks  using  a:  (b) 
Matched  filter,  (c)  Phase-only  filter,  (d)  lEJTC  at  m  =  10'^, 
(e)  lEJTC  at  m=10\  (f)  lEJTC  at  m=10^ 

(4)  (a)  The  input  plane  with  tanks  used  for  the  simulation 

Correlation  peaks  intensity  plots  for  the  two  tanks  using  a:  (b) 
Matched  filter,  (c)  Phase-only  filter,  (d)  lEJTC  at  m  =  10*^, 
(e)  lEJTC  at  m  =  10^  (f)  lEJTC  at  m  =  10^ 
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Chapter  3 


Speckle  Velocimeter 


3.1.0  Abstract 

In  this  chapter  we  propose  a  speckle  velocimeter  which  is 
based  on  using  a  real-time  hologram  in  four  wave  mixing  as  a  time- 
integrative  correlator.  The  theory  of  the  speckle  velocimeter  has 
been  developed  for  the  time  correlation  between  the  far-field 
spectrum  of  light  scattered  from  the  diffuser  and  reference  wave 
which  is  doppler  shifted.  Our  theoretical  derivation  shows  that  the 
output  from  the  time-integrative  correlator  is  proportional  to  the 
square  of  the  velocity.  The  proportionality  factor  is  dependent  on 
three  factors:  (a)  a  complicated  function  of  the  ratio  between  the 
surface  roughness  and  the  granularity;  (b)  the  input  beam  ratios; 
and  (c)  the  coupling  coefficient  of  the  photorefractive  media. 
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3.2.1  INTRODUCTION 

Speckle  metrology  (1)  and  velocimetry (2-4)  have  been  applied 
to  numerous  applications  in  solid  and  fluid  mechanics.  These 
applications  range  from  fluid-flow  research,  wind  tunnel  testing, 
flow  control  system  and  combustion.  There  is  substantial 
j_-terature  describing  various  designs,  theory  and  operational 
mechanisms  of  various  speckle  velocity  meters (2-8)  .  However,  the 
operational  mechanism  of  these  devices  can  be  divided  into  two  main 
categories.  The  first  category  is  the  interferometric  (9,10,11)  and 
the  second  is  time  correlation  (6,7,8).  The  interferometric 
techniques  have  been  applied  by  using  dovible  exposure  techniques  in 
either  a  real-time  hologram  or  a  photographic (9, 10, 11)  film  or  by 
using  an  array  of  CCD  (13)  detectors. 

The  implementation  using  photographic  films  is  not  real-time 
because  a  chemical  process  is  required  to  develop  the  hologram. 
However,  the  implementation  using  real-time  holography  or  a  CCD 
array  is  also  not  real-time  because  it  still  requires  substantial 
digital  processing  which  compromises  the  advantages  sought  by  the 
optical  signal  processing. 

The  other  approach  which  measures  the  autocorrelation  function 
assumes  that  when  a  coherent  light  illuminates  a  moving  rough 
surface  the  generated  speckle  pattern  fluctuates  accordingly. 
Therefore,  if  a  photodetector  records  the  variation  of  the  light 
intensity,  that  should  be  a  measure  of  the  velocity  of  the 
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diffuser.  The  simplest  method  to  determine  the  velocity  is  to 
measure  the  autocorrelation  function  of  the  speckle  pattern. 

The  measurement  of  the  autocorrelation  function  has  been  used  to 
determine  the  speed  of  linearly  moving  and  rotating  objects.  The 
drawback  is  that,  these  systems  measure  the  fluctuation  only  at  one 
point  in  space  which  does  not  give  enough  sampling  points.  The 
second  drawback  is  that  these  systems  are  is  not  immune  to 
background  vibrational  noise.  To  overcome  this  problem,  Moran,  et 
al.,  described  a  new  way  to  measure  the  correlation  function  using 
a  phase-locked  speckle  electronic  pattern  (6,7). 

The  new  system  used  an  array  of  CCD  detectors  in  order  to 
measure  the  fluctuation  over  the  whole  field.  The  background  noise 
was  isolated  by  using  phase-lock  homodyne  or  heterodyne  detection 
,  in  order  to  demodulate  the  fluctuation  in  the  speckle  pattern 
from  the  background  noise.  The  reference  signal  was  generated 
electronically . 

This  system  has  demonstrated  the  ability  to  perform  with  high 
sensitivity  even  the  presence  of  vibrational  background  noise  whose 
level  would  destroy  the  performance  of  conventional  non-isolated 
speckle  velocimeters . 

Recently,  we  demonstrated  the  first  optical  lock-in  amplifier 
which  operates  on  the  principle  of  time  correlation  using  four-wave 
mixing  in  real-time  holography (13 , 14 , 15)  .  This  optical  lock-in 
amplifier  has  been  used  by  Lockheed  for  analysis  of  vibration  (16)  . 
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In  manner  similar  to  Moran,  et  al.  our  lock-in  amplifier  can  be 
used  as  a  speckle  velocimeter.  The  advantage  here  is  that  the  new 
design  is  simpler  to  implement  and  more  cost  effective  by  using  a 
large  plate  of  photorefractive  polymers  (17).  Further  more  our 
proposed  device  gives  a  direct  measurement  of  the  speed  without  any 
further  processing. 

In  this  chapter  we  will  present  (a)  the  theory  of  light 
scattered  out  of  a  diffuser  moving  with  longitudinal  motion,  (b)  the 
description  of  the  device,  (c)  the  theory  of  the  two-dimensional 
time  correlation  detector  using  photorefractive  media  and  (d)  the 
time  correlation  between  speckle  from  a  diffuser  moving  with 
longitudinal  motion  and  doppler-shifted  light. 
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3.2.2  SPECTRUM  OP  LIGHT  SCATTERED  FROM  A  LONGITUDINALLY  MOVING 


DIFFUSER 

Various  theories  in  the  past  have  been  developed  to  describe 
the  far  field  spectrum  of  light  scattered  from  a  diffuser.  The 
theories  have  already  been  developed  for  transverse  (18) , 
longitudinal (19)  and  rotational  motions  of  diffusers  (20,21  ). 

Since  we  are  going  to  develop  the  theory  of  speckle  velocimeter  of 
diffuser  with  a  transverse  motion,  therefore  we  are  going  to 
represent  the  theory  of  light  scattered  from  a  moving  diffuser  with 
transverse  motion  (18)  . 

The  far  field  spectrum  S“  of  light  scattered  from  a 
longitudinally  moving  diffuser  of  aperture  A  is  given  by: 

ru')  =—  (-^^)  2Q'"c:(Q^',iOS',(Q") 

2it  cr  ^  ^ 


where  Cj,  is  the  two-dimensional  spatial  Fourier  transform  of  the 
diffuser  correlation  function.  5,(0'')  is  the  spectrum  of  the 
incident  light  on  the  diffuser,  r  is  the  distance  between  the 
aperture  and  the  observation  point, c  isthe  speed  of  light,  6  is  the 
angle  that  the  direction  u'  makes  with  the  z  axis,  u'per  is  the 
projection  ,  considered  to  be  a  two-dimensional  vector  unit  u'on 
the  aperture  A.  n',n"  are  the  frequencies  of  the  light  scattered  by 
and  incident  on  the  diffuser,  respectively.  D'and  D”  are  related  by 
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the  following  formula: 


Q/'=Q/(l-pu/)  =Q/(i-pcos<i>) 

in  which  jS=v/c,  and  <p  is  the  angle  between  u  and  v.  See  Fig  1(a) 
and  1(b)  for  further  illustrations 

The  Fourier  transform  of  the  correlation  function  is  given  by: 
^(Q,f)  =-^j*C^(Q,Q')  expi-if.q') 


where  p  is  the  vectorial  difference  between  two  arbitrary  points  in 
the  diffuser,  while  the  diffuser  correlation  function  is  given  by 

Cc(Q  :e)  *exp(T)  (fl)  (l-g(Q)) 


g  is  the  normalized  correlation  function  of  the  diffuser  height  and 
is  defined  as: 


is  the  variance  of  the  plate  thickness  and  is  given  as: 
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<j|=<  IhiQ)  -<i2(Q)  >]  2 


For  a  diffuser  which  obeys  the  Poisson  distribution (22)  (i.e.  in 

a  diffuser  of  width  W  and  cell  size  l^j,  the  probability  of  having 

T  cells  within  the  interval  W  is  given  by: 

Pj.=  (P//I^)’’exp(-JV/i^)  r=0,l,2 . 

0  otherwise 

The  normalized  correlation  function  for  a  binary  or  multilevel 
diffuser  with  a  Poisson  distribution  is  given  by: 

g{Q)=l--^  when 

0  when  Q,>lc 


For  a  ground  glass  diffuser  with  Poisson  distribution,  the 
normalized  correlation  function  is  given  by: 


g(Q)  =exp-8/I„=l— 

i|  I'  ll 


Substituting  Eq.  (8)  in  Eq.  (4)  yields: 

C^(Q,  q)  =exp  (n(Q)  -1)  -^e]  2= 
exp-^[(n{Q)-l)-^]2 


29 


For  a  ground  glass  diffuser  the  additional  terms  in  Eg.  (9)  can  be 
included.  By  calculation  of  the  Fourier  transform  of  Eg.  3  and 
substitution  into  Eg. (1)  yields  the  following  result: 


N 


n(Q")  -1 


S^iQ")  e3g)- 


tl,sine]2 


2  [a^(n(Q//)-l)  (l-pu/)]2 


where  N  is  a  freguency- independent  normalization  factor,  defined 
as: 


AlpCos^Q 

N= - - - - 

[2nro^(l-puj')  ]  ^ 


The  interesting  feature  of  Eg.  (11)  is  that  the  Doppler  shift  of  the 
scattered  light  is  mostly  dependent  on  v/c  cos  <f),  where  0  is  the 
angle  between  uper  and  v  .  The  weighting  factor  of  the  Doppler 
shift  (or  the  power  spectrum)  at  a  certain  angle  6  is  dependent 
mainly  on  the  ratio  between  1^  and  cr,,.  The  dependence  in  the 
previous  form  looks  to  be  very  logical  for  the  following  reason. 
When  light  is  scattered  out  of  a  diffuser,  the  light  can  be 
scattred  into  the  far  field  in  any  angle  between  6=0  up  to  *6=90. 
According  to  the  doppler  shift  theory,  the  freguency  of  the 
scattered  light  at  a  certain  position  should  be  dependent  only  on 
the  speed  of  the  source  relative  to  the  point  of  observation  (i.e. 
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on  6 ,  0  and  v)  .  The  dependence  of  the  power  spectra  on  l^.  and 
because  the  light  scattering  dependent  only  in  the  diffuser  cell 
characteristic  even  further  on  the  ratio  between  l^j  and  This 
because  a  larger  ah  will  make  the  variation  in  the  phase  of  the 
scattered  light  faster,  while  a  larger  Ic  deminishaes  the  variation 
in  the  phase  of  the  scattered  light. 
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3.2.3  THE  TIME  INTEGRATIVE  CORRELATOR 

Consider  three  beams  of  the  same  wavelength  A,,  A^,  and 

intersecting  within  a  photorefractive  medium,  as  shown  in  Fig  (2)  . 
A,  and  A4  are  coherent  so  they  combine  to  induce  a  hologram  within 
the  material;  while  Ag  are  coherent  within  themselves,  but  not  with 
A^  and  A4.  Mathematically  the  amplitudes  are 

Aj^=E^{t)  exp  [i  ik^r-a^t)  ] 

^2=^2  ( t)  exp  [i  (icjr-coa  t)  ] 

A^=E^  ( t)  exp  [i  t)  ] 

In  the  steady  state  case  (i.e  beams  A^  and  A^  not  varying  in  time 
)  a  hologram  is  written  and  a  phase  conjugate  is  generated.  This 
phase  conjugation  is  counter  propagating  with  respect  to  beam  A^. 
Treatment  of  the  time  varying  input  necessitates  the  solution  of 
the  photorefractive  space  charge  field  in  the  dynamic  regime.  For 
the  sake  of  simplicity  we  will  assume  that  the  total  intensity  of 
the  beam  incident  on  the  cirystal  is  constant  (i.e  any  amplitude 
variation  in  the  input  beams  is  considered  small  compared  to  the 
summation  of  the  intensities  of  all  the  beams) .  The  development  of 
the  space  charge  field  (  within  the  crystal  is  described  by 

^5C  -  Y 

dt  X  X  Iq 
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where  r  is  the  response  time  of  the  grating  formation  ,  y  is  the 
coupling  coefficient  of  the  material,  and  I,  is  the  total  intensity 
incident  on  the  crystal.  Equation  (14)  assumes  that  the  readout 
beam  is  sufficiently  weak  so  that  it  does  not  effect  the  grating. 
The  general  solution  of  Eq.  (14)  is 


_c 


Jo  X 


Equation  (15)  can  be  treated  by  assuming  that  A^(t)  and  A4(t)  are 
signals  with  Fourier  transforms  such  that 


JE 


i  pY^i(t)A4  (t) 

Jo  X 


Equation  (16)  can  be  treated  by  assuming  that  A^(t)  and  A4(t)  are 
signals  with  Fourier  transforms  such  that 


A^(C)  exp(-io)'t)  do' 

oo 

A^{t)=js^  (o")  exp  (-io"t)  do" 


Substituting  Eq.  (17)  into  Eq.  (16)  yields 


33 


-  (q'-co'O  +— 

X 

exgl-i  (w'-w^O  t]  dtii'dai" 


In  the  diffusion  limit  t  is  a  real  number.  With  a  dc  electric  field 
applied  to  the  crystal,  t  becomes  complex: 


— +io). 


This  yields 


5,  (0)0  5^  (o>^0 
(G)'-0)"+Ci)P  + 


exp  t-i  ((i)'-(i)"+«y)  -|]  dw'  din" 


3.2.4  TIME  CORRELATION  WITH  SPECKLE 


using 


Figure  3  shows  a  schematic  diagram  of  the  speckle  velocimeter 
a  photorefractive  medium.  A  laser  beam  incident  on  the 
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acousto-optic  modulator  Is  divided  into  two  beams:  and  A^.  Beam 
A4  is  the  deflected  doppler-shifted  component  of  the  beam,  while 
beam  A,  is  the  transmitted  component.  Beam  A^  passes  through  a 
moving  diffuser  thereby  generating  a  nonstationary  speckle  pattern 
with  doppler  shift. 

The  light  from  the  two  beams  is  combined  together  in  a 
photorefractive  crystal.  The  recombined  light  can  write  a  grating 
as  long  as  the  doppler  shift  in  beam  A,  is  near  the  doppler  shift 
of  some  of  the  light  scattered  by  the  diffuser.  Beam  Aj  is  a  read¬ 
out  beam  of  the  grating.  The  scattered  light  of  the  grating  is 
directed  by  the  beam  splitter  into  the  detector,  which  can  also  be 
a  CCD  array.  The  measured  light  is  analyzed  to  detect  the  velocity 
of  the  diffuser. 

Let  us  consider  that  the  reference  beam  is  given  by: 


A,  ( t)  =A,e 

or,  in  the  spectrum  domain  it  is  given  by 


For  the  sake  of  simplicity,  it  is  possible  to  assume  that  there  is 
no  dispersion  in  the  index  of  refraction  of  the  diffuser. 
Accordingly,  let  us  assume  that  r\(n")=Tig  and  that  the  incident 
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light  into  the  diffuser  is  monochromatic  at  the  frequency  n”, 
i.e. , 

iui")  =5 


By  replacing  n”  by  o"-n''  and  substituting  Eqs  (11)  ,  Eq  (23)  and 
(24)  in  to  Eq.  (2),  it  is  possible  to  show  that. 


_  A^A,N  (i^sine)^ 

[22^-1]^  2((o^(n^-l)  (1-Pu^)]2 


] 


x[-^(  — +io)J  ] 


-j  +  — 


Let  us  define  a  new  parameter,  which  we  shall  call  the  diffuser 
cell  characteristic  parameter  (DCCP)  ,  a^.  This  parameter  is  defined 
as: 

Substituting  (DCCP)  in  Eq.  (25)  yields  the  following  format: 


The  integrated  intensity  of  the  output  is  given  by: 
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This  integral  is  very  difficult  to  evaluate  analytically.  However, 
it  is  still  possible  to  do  some  approximations  in  order  to  evaluate 
it.  We  can  consider  that  the  second  term  in  Eq.  (27)  is  dependent 
on  6,  because  the  value  of  v/c  is  always  very  small,  While  the 
third  term  in  the  same  equation  can  be  expanded  into  a  Taylor 
series  because  in  the  lock-in  approximation: 

1  (tOj^-Q^+cOg.)  xl<l 


For  evaluating  the  whole  output,  it  is  necessary  to  evaluate  the 
integration  of  the  squares  of  the  second  and  third  terms  of  Eq. 
(27)  . 

By  neglecting  the  dependence  on  i  in  the  square  of  the  integral  of 
Eq.  (27),  then  this  integral  can  be  approximated  as: 


While  the  integration  of  the  square  of  the  third  term  for  the  first 
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two  orders  of  the  Taylor  expansion  is  given  by: 
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Let  us  assiame  that  the  frequency  of  the  light  scattered  off  the 
acousto-optic  transducer  is  given  by: 

where  v,  is  the  speed  of  the  propagation  of  the  sound  in  the 
acousto-optic  modulator,  is  the  frequency  of  the  laser  source 
which  is  also  equal  to  n-  (i.e  o,=n«) .  The  integration  of  first 
term  in  the  series  in  Eq.  (32)  is  equal  to: 

[l+(<0j^x)2]  [n(a)^P)2+2x(w^Pi-C0y)2] 

The  integration  of  the  second  term  in  the  series  of  Eq.  (32), 
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includes  terms  to  the  power  four  in  )31  and  ^  and  cross  terms  to 
the  power  two  in  ppl,  therefore,  these  latter  terms  are  neglible 
compared  to  the  first  term. 

According  to  the  integral  evaluation,  the  total  intensity  at  the 
output  is  given  by  the  following  formula: 

^cut=  ^ )  5^  ( [it(w^p2) +27:(o)„Pi-o)p2] 


where  r  is  the  coupling  coefficient  in  the  intensity.  It  is  evident 
from  this  equation  that  the  integrated  intensity  of  the  output  is 
dependent  on  the  square  of  the  speed.  This  fact  makes  the  current 
device  very  easy  to  implement  experimentally  in  the  sense  that 
there  is  no  need  for  extra  processing  of  the  output.  A  simple 
electronic  circuit  which  takes  the  square  root  of  the  current 
should  gave  a  direct  measurement  of  the  speed. 

So  far,  we  have  developed  the  theory  of  speckle  velocimetry 
based  on  the  far  field  spectrum  from  a  moving  diffuser .However,  the 
current  theory  should  be  valid  for  other  applications  in  fluid 
mechanics  and  biomedical  engineering.  The  only  difference  will  be 
in  the  term  H(ajj)  .  This  factor  should  be  more  characteristic  for 
the  speckle  pattern  which  is  scattered  out  of  the  diffuser,  and 
would  have  to  be  recalculated  for  each  specific  problem. 
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3.2.5  Experimental  Results. 

In  the  past  we  performed  some  initial  experiments  which 
justifies  the  concept  of  speckle  velocimetry  using  the 
photorefractive  time  integrative  correlator.  Fig.  4  shows  the 
experimental  arrangement  which  was  used.  This  is  essentially  a 
standard  configuration  for  four-wave  mixing  modified  by  two  PZTs. 
The  two  PZTs  are  driven  by  two  separate  function  generators .  The 
source  was  a  HeNe  laser  operating  in  a  multilongitudinal  mode  with 
an  output  power  of  30  mW  and  wavelength  1=632. 8nm.  The  beam  was 
expanded  to  a  diameter  of  2mm.  Beam  splitter  BS,  divided  the  source 
beam  into  two  beams:  reference  beam  and  the  signal  beam. 

The  readout  beam  A2  was  generated  by  retroreflecting  A^  from  mirror 
M3.  A  lens  (f=100mm)  was  used  in  beam  4  to  image  the  input  object 
plane  to  the  BSO.  The  measured  intensities  of  beams  I,,  12,14  were 
10  mW,  5mW,  and  8mW,  respectively.  In  this  experiment  we  used  a 
square  wave  modulation  instead  of  using  ramp  phase  modulation.  We 
chose  this  kind  of  modulation  in  order  to  maximize  the  erasing 
tendency  of  the  modulation  in  a  coherent  erasure  technique.  (This 
may  change  the  results  slightly,  but  not  significantly,  because 
higher  orders  of  locking  will  not  contribute  that  much  to  the  total 
intensity  at  the  output.  For  static  objects  under  continuous 
flipping  of  the  phase,  the  phase  conjugate  beam  intensity  was 
averaged  to  zero.  However,  for  moving  diffusers  the  phase  conjugate 
was  not  averaged  to  zero;  reappearing  under  the  constraints  of  time 
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correlation  detection.  Figure  5(a)  is  the  conjugate  image  of  the 
diffuser.  Fig  5(b)  shows  the  results  under  of  the  erasure  by  a 
square  phase  modulation.  When  the  diffuser  translated  in  any 
direction,  a  phase  conjugate  image  appeared  as  shown  in  Fig  6c.  The 
irradiance  of  the  image  is  very  dependent  on  the  speed  of  the 
motion.  Slow  motion  was  very  faint,  with  increasing  speed  the 
output  irradiance  increased  until  the  frequencies  of  the  variation 
in  the  input  image  no  longer  matched  the  reference  beam 
frequencies.  These  qualitative  results  prove  the  operational 
principle  of  speckle  velocimetry.  However,  further  quantitative 
work  on  this  aspect  needs  to  be  done. 

3.2.6  Conclusion: 

We  have  presented  the  theory  of  speckle  velocimetry.  Our  theory  is 
based  on  the  use  of  the  photorefractive  time-integrative  correlator 
and  diffuser  which  is  moving  with  a  constant  speed  .  In  our  theory 
we  found  that  the  output  intensity  is  proportional  to  the  square  of 
diffuser  velocity.  In  contrast  to  many  previously  proposed  designs, 
the  velocity  is  easily  extracted  without  the  necessity  for  any 
further  processing  of  the  output.  The  new  current  theory  should 
open  the  possibility  of  considering  the  use  of  this  device  in  fluid 
mechanics  problems  and  biomedical  applications. 

A  priliminary  experimental  result  which  was  obtained 
previously  justifies  the  current  theory.  However,  further 
experimental  work  needs  to  be  done  to  understand  the  theory  for 
other  cases. 
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3.2.7  Figure  Captions. 

Fig.  1  (a)  Illustrative  scheme  to  determine  the  far  field  spectrum 
of  light  scattered  from  a  diffuser,  (b)  Illustrative  scheme  to  show 
the  projection  of  U  in  the  plane  of  motion  and  the  angle  between 
uper  and  v. 

Fig.  2  A  schematic  diagram  to  illustrate  the  time  integrative 
correlator  with  piezo-mirrors. 

Fig  3.  (a)  A  diagram  of  a  speckle  velocimeter  using  an  acousto¬ 
optic  modulator  which  serves  as  a  beam  splitter  as  well  as  a 
doppler  shift  device,  (b)  The  geometrical  arrangement  to  produce 
the  effect  of  scattered  light  from  transverse  motion  of  diffuser 
from  a  rotating  diffuser. 


Fig  4.  A  schematic  diagram  of  a  speckle  velocimeter  in  the  imaging 
geometry . 

Fig  5.  Illustration  of  speckle  velocimetry  with  transverse  motion 
of  the  diffuser. (a)  The  phase  conjugate  image  with  a  static 
diffuser,  (b)  With  square  phase  modulation  operate  on  PZTl.  (c)  The 
regenerated  phase  conjugate  with  transverse  motion  of  the  diffuser. 
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chapter  4 


Homodyne  and  Heterodyne  Imaging  in  a 
Scattering  Medium 


4.3.0  Abstract 

In  this  paper  we  introduce  two  dimensional  homodyne  and 
heterodyne  techniques  for  imaging  objects  in  the  presence  of  a 
scattering  medium.  These  techniques  are  categorized  into  two  forms 
according  to  their  operational  principles.  The  first  is  based  on 
the  Doppler  shift  difference  in  light  scattered  from  objects  of  two 
different  textures  or  an  opaque  object  and  a  textured  scattering 
medium.  The  second  technique  is  based  on  the  principle  of  photon 
density  waves,  but  incorporates  2-d  holographic  detection.  In  both 
cases  a  two-dimensional  real-time  holographic  phase-sensitive 
detector  is  proposed  to  pull  the  images  from  noise.  We  report 
the  initial  demonstration  of  pulling  signals  out  of  noise  based  on 
the  Doppler  shift  difference.  Enhancements  of  signal-to-noise  in 
the  order  of  50  have  been  achieved. 
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4.3.1  INTRODUCTION: 


Imaging  through  a  scattering  medium  or  imaging  an  object 
within  a  scattering  medium  is  one  of  the  most  challenging  problems 
in  optical  signal  processing.  It  has  significant  importance  in 
tomography  and  image  formation  in  biological  tissues.  Various 
methods  have  been  utilized  to  solve  this.  However,  one  of  the  most 
popular  approaches  is  based  on  the  principle  of  the  first  arriving 
light.  In  this  technique,  a  time  gate  is  used  to  separate  the 
information-carrying  light  (the  first  light)  from  the  noisy  light 
(the  subsequent  light) .  The  first  realization  of  this  scheme  was 
demonstrated  by  Mantic  and  Duguay(l) .  In  their  specific 
implementation  the  hologram  was  used  as  a  time  gating  processor. 
Since  then  ,  different  time  gating  techniques  have  been  used,  for 
example,  electronic  holograms(2)  and  real-time  holograms (3) .  Time 
gating  has  also  been  implemented  using  Kerr (4)  and  Raman  cells (5) . 
Further  improvement  in  terms  of  signal-to-noise  ratio  have  been 
achieved  by  adding  spatial  filtering  techniques (6) . 

One  drawback  in  the  first  light  techniques  is  they  require 
ultrafast  pulses  and  sophisticated  instrumentation  for  fast 
imaging.  These  facts  make  the  above  approaches  expensiye  to 
implement.  Therefore,  we  propose  an  alternative  technique  which 
uses  a  CW  laser  associated  with  homodyne  and  heterodyne  detection. 
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2.3.2  Background: 

In  this  chapter  we  introduce  two-dimensional  homodyne  and 
heterodyne  techniques  for  imaging  objects  in  the  presence  of  a 
scattering  medium.  Imaging  in  scattering  media  is  important  for 
biomedical  applications  (for  example,  imaging  tamors  within  a 
healthy  tissue)  or  transferring  images  within  atmospheric 
turbulence.  Our  techniques  are  categorized  into  two  forms  according 
to  their  operational  principle.  The  first  is  based  on  the  Doppler 
shift  difference  in  light  scattered  from  objects  of  two  different 
textures  or  an  opaque  object  and  a  textured  scattering  medium.  We 
also  show  that  the  basic  principle  of  the  Doppler  shift  difference 
technique  and  of  first  light  are  similar.  The  advantage  is  that  the 
Doppler  technique  does  not  require  a  high  power  pulsed  laser. 
Rather,  CW  lasers  with  modulation  in  the  range  of  MHz  to  GHz  are 
used.  The  second  technique  is  based  on  the  principle  of  photon 
/density  waves (7 ) ,  but  incorporates  2-d  holographic  detection.  In 
both  cases,  a  two-dimensional  real  time  holographic  phase  sensitive 
detector  is  used  to  pull  the  images  from  noise  (8,9,10,11,12). 
Opto-electronic  devices,  based  on  smart  pixelated  arrays,  can  also 
be  used  for  this  purpose. 

In  the  first  light  technique  a  narrow  pulse  of  light  is  passed 
through  a  scattering  medium.  The  first  emerging  light  is.  less 
severely  scattered  and  therefore  forms  the  best  image,  while  the 
subsequent  light  is  considered  noise.  The  problem  then  reduces  to 
extracting  signal  from  noise  by  time  gating. 
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In  the  Doppler  difference  technique,  cw  light  is  scattered  by 
moving  media  of  different  textures  (having  different  sized 
scattering  centers) .  It  is  well  known  in  statistical  optics,  that 
scattered  light  from  a  random  media  moving  with  a  constant  velocity 
will  cause  Doppler  shifting ( 13 , 14) .  This  is  true  if  the  motion  is 
along,  or  perpendicular  to,  the  direction  of  propagation  of  the 
illuminating  light.  An  object  with  large  scattering  centers 
produces  a  relatively  small  angle  of  scatter  and  hence  generates  a 
small  Doppler  shift.  In  contrast,  small  scattering  centers 
generate  larger  angles  of  scatter  and  therefore  larger  Doppler 
shifts.  For  example,  consider  an  object  of  one  texture  embedded 
inside  a  medium  of  another  texture,  as  in  the  case  of  a  tumor  in 
healthy  tissue.  Each  texture  will  contribute  different  Doppler 
shifts.  Clearly,  homodyne  and  heterodyne  techniques  will  be  useful 
in  separating  the  light  scattered  from  each  of  these  two  textures. 

Large  motions  of  an  object  are  not  very  practical.  For 
example,  if  a  tumor  and  healthy  tissue  are  moving  simultaneously, 
the  image  would  be  smeared.  In  the  case  of  a  massive  object,  this 
problem  can  be  overcome  by  vibrating  only  the  scattering  medium, 
for  example,  by  launching  an  acoustic  wave.  For  non-sinusoidal 
vibrational  motion  of  a  random  scattering  medium,  the  Doppler  shift 
spectrum  shows  multiple  peaks  at  harmonics  of  the  fundamental 
vibrational  frequency,  each  with  its  own  characteristic  bandwidth 
depending  on  the  statistics  of  the  scattering  medium.  Objects  with 
coarser  texture  will  of  course  scatter  light  with  less  bandwidth. 
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This  light  can  be  separated  from  the  other  light  using  frequency 
selective  detection. 

In  previous  publications,  we  reported  on  a  two-dimensional 
real-time  holographic  optical  lock-in  detector(8,9,10) .  It  was 
also  demonstrated  that  it  is  possible  to  detect  moving  texture.  As 
in  all  heterodyne  detection  schemes,  maximxim  detection  of  moving 
texture  occurs  when  the  center  of  the  Doppler-shifted  frequency  of 
the  scattered  light  from  the  random  medium  is  equal  to  the 
modulation  frequency  of  the  reference  optical  signal (15, 16, 17) . 
Nearby  frequencies  are  also  recovered  and  contain  additional 
information  as  in  all-optical  baseband  demodulation(18,19) . 

4.3.2  DEVICE  DESCRIPTION: 

Figure  1  shows  a  schematic  diagram  of  a  homodyne  and 
heterodyne  scheme  for  imaging  within  the  scattering  medium  (A)  for 
an  object  behind  the  scattering  medium  (B)  for  object  imbedded  in 
the  scattering  medixim  and  (C)  for  an  embedded  object,  but  with  the 
light  itself  modulated,  rather  than  the  scattering  medium.  Cases 
(A)  and  (B)  are  based  on  the  principle  of  the  difference  in  Doppler 
shifting  while  (C)is  based  on  the  principle  of  photon  density 
waves,  and  can  be  extended  to  the  case  of  multiple  light  sources 
with  fixed  phase  shifts. 

In  all  cases  a  two-dimensional  lock-in  detector  is  used.  This 
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two-dimensional  lock— in  detector  can  be  a  real-time  hologram 
operating  with  various  configurations  such  as  four-wave 
mixing (8, 9 , 10) ,  two-wave  mixing(ll)  or  a  two-input  version  of  a 
self— pumped  phase  conjugator (12) .  In  all  of  these  configurations, 
one  of  the  input  beams  has  the  signal  and  the  other  beam  has  the 
modulated  reference.  Another  way  to  implement  the  two-dimensional 
lock-in  detector  is  to  use  a  smart  pixelated  array  where  the 
homodyne  and  heterodyne  detection  can  be  controlled  by  an  internal 
clock. 

In  (A)  the  object  is  set  behind  the  scattering  medium.  In  this 
case  laser  light  passes  first  through  the  object  and  then  through 
the  scattering  medium.  An  oscillation  is  generated  by  an  acoustic 
vibrator  attached  to  the  scattering  medium  or  some  other  acoustical 
transducer.  In  this  case  the  zero  Doppler-shifted  light  corresponds 
to  the  light  which  carries  the  information  on  the  signal,  while  the 
Doppler  shifted  light  corresponds  to  the  scattered  light. 
Therefore,  homodyne  detection  for  the  non-shifted  light  should 
recover  the  signal.  Homodyne  detection  on  the  non-Doppler-shifted 
light  does  not  require  a  modulated  reference.  Thus,  a  simple  CCD 
array  or  a  time  averaging  hologram  is  sufficient  for  recovery  of 
the  signal  from  the  noise.  Case  (B)  is  similar  to  case  (A)  except 
that  the  object  is  embedded  in  the  scattering  medium.  Simultaneous 
vibration  of  an  object  and  a  scattering  medium  of  different 
textures  will  cause  each  to  scatter  light  with  a  different  Doppler 
shift.  Therefore,  using  a  homodyne  detection  technique  and  proper 
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choice  of  reference  frequency,  it  is  possible  to  lock  in  on  the 
maximum  of  the  frequency  spectnim  of  either  the  signal  or  the 
noise,  with  some  bandwidth  tolerance.  (17,18,19) 


The  last  configuration  for  detecting  a  signal  in  a  noisy 
environment  is  based  on  the  photon  density  wave  experiments.  In 
this  case,  the  demodulated  light  carries  the  phase  information  of 
the  photon  density  wave  after  it  has  diffracted  off  the  object  and 
propagated  through  the  scattering  medium.  Implementation  of  this 
scheme  with  a  single  (scanning)  photodiode  (1-d)  has  been  described 
by  others  [7].  The  improvement  described  here  is  to  perform  2-'d 
parallel  optical  demodulation  of  the  scattered  light.  The  basic 
concept  of  photon  density  wave  imaging  relies  on  phase  sensitive 
detection  of  an  amplitude  modulated  laser  diode.  We  have  shown 
:  that  the  phase  of  the  modulation  of  a  laser  diode  beam  can  be 
directly  measured  using  all-optical  phase  sensitive  detection  for 
modulation  frequencies  up  to  several  GHz.  This  phase  measurement 
was  accomplished  by  varying  the  path  length  transverse  by  the 
modulated  light  and  monitoring  the  optical  lock-in  demodulated 
signal  [20] . 

4.2.4  EXPERIMENT: 


Two  experiments  have  been  performed  to  demonstrate  the 
principle  of  imaging  an  object  behind  the  scattering  medium.  In  the 
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first  experiment  a  CCD  camera  is  used,  and  in  the  second,  a  four 
wave  mixer  with  a  photorefractive  BaTi03  crystal.  Both  schemes  are 
shown  in  Fig  2(A)  and  2(B)  respectively  .  As  shown,  in  both  cases, 
light  reflected  of  the  resolution  chart  (the  object  )  is 
transmitted  through  nozero  mean  diffuser  which  is  mounted  on  an 
acoustic  vibrator,  driven  by  sawtooth  voltage.  The  diffuser 
consists  of  a  spray-painted  microscope  slide.  A  lens  focal  length 
30cm  was  used  to  image  the  resolution  chart  in  to  the  camera  (fig 
2A)  .  Similarly  in  fig.  2B,  the  CCD  camera  was  replaced  by  four  wave 
mixer  with  BaTi03  crystal,  lens  Ij  was  used  to  image  the  phase 
conjugate  output  into  the  CCD  camera.  The  write  beams  (1  and  4) 
were  extraordinary  polarized,  and  the  read  out  beam  was  ordinary 
polarized.  These  polarizations  were  chosen  to  prevent  beam  fanning 
of  the  crystal  during  the  read  out  process,  and  to  simplify 
separation  of  the  phase  conjugate  and  the  light  scattered  by  the 
crystal  surface  through  using  a  polarizer. 

Figure  3  shows  the  experimental  results  for  imaging  within  a 
scattering  medium:  (A)  Is  the  orginal  image  (B)  the  light  incident 
on  the  CCD  after  transmission  through  the  image  and  a  non-moving 
diffusser,  and  (C)  the  light  when  the  diffuser  is  vibrated  at  a 
frequency  of  lOHz.  The  improvement  in  signal  to  noise  ratio  was  20. 
In  Fig  3(B)  and  3(C),  the  experiment  is  repeated  using  the 
photorefractive  four-wave  mixer  and  a  vibration  frequency  of  .... 
This  frequency  is  less  than  the  previous  one  because  of  the  long 
time  integration  of  the  BaTio3 .  In  this  case,  the  SNR  is  improved 
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from  50  (Fig  3B)  to  (Fig  3D). 


The  results  showed  better  noise  removal  using  the 
photorefractive  crystal.  This  can  be  attributed  to  two  reasons.  The 
first  reason  is  due  to  amplification  of  the  signal  beam  by  the 
crystal.  The  second  is  due  to  rejection  of  the  scattered  light  by 
the  crystal,  compared  to  simply  averaging  out  the  spatial  content 
of  the  signal,  as  is  the  case  with  the  camera  only. 

4.3.5  CONCLUSION: 

In  conclusion,  we  have  proposed  homodyne  and  heterodyne 
techniques  for  imaging  in  a  scattering  medium.  The  first  technique 
is  based  on  the  principle  of  the  Doppler  shift  difference  of  the 
light  which  is  scattered  from  an  object  vs  light  scattered  by  the 
surroundings.  We  showed  that  the  principle  of  Doppler  shift 
difference,  and  the  principle  of  first  arriving  light  are  similar. 
However,  in  contrast  to  first  light,  which  requires  ultrafast 
pulses,  here  modulated  CW  can  be  used.  We  also  experimentally 
demonstrated  the  basic  principle  for  a  simple  case  of  an  object 
outside  of  the  scattering  medium.  The  second  is  an  extension  of  the 
principle  of  photon  density  waves,  using  a  two-dimensional  (non¬ 
scanning)  detection  scheme. 

Many  nonlinear  media  can  be  used  for  the  processes  described. 
The  preference  is  for  high  gain,  low  noise  materials.  The 
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materials  which  we  recommend  is  ferroelectric  photorefractive 
materials  such  as  BaTiOj  with  Rh  doping  [21]  in  self-bending 
geometry  [22].  This  material  has  the  largest  coupling  coefficient 
observed  so  far  and  works  in  the  wavelength  range  of  primary 
interest  for  biological  tissues.  Other  materials  such  as  thin  film 
polymers  [23]  can  be  used  with  multistage  amplification  because  one 
thin  film  is  not  enough.  Resonant  systems,  such  as  atomic  vapors 
[24],  can  also  be  used  where  high  speed  response  at  low  light 
levels  are  required.  In  all  of  these  geometries  we  can  use  a 
diverging  pump  beam  to  ensure  high  resolution  gain  [25]. 
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4.3.6  FIGURE  CAPTIONS: 


Fig:  1:  A  schematic  diagram  which  illustrates  the  principle  of 
imaging  in  a  scattering  medium  based  on  using  two-dimensional 
phase  sensitive  detector  Figures  (A)  and  (B)  are  based  on  the 
principle  of  the  Doppler  shift  difference:  (A)  for  object  impeded 
behind  the  scattering  medium;  (B)  For  an  object  embedded  in  the 
scattering  medium;  (C)  Is  based  on  the  principle  of  photon  density 
waves  for  an  object  which  embedded  either  behind  or  within  the 
scattering  medium. 

Fig  2:  The  experimental  set-up  used  for  imaging  an  object  imbedded 
behind  a  scattering  medium:  (a)  is  based  on  using  the  CCD  camera  as 
the  two-dimensional  phase  sensitive  detector  (B)  base  on  using  a 
photorefractive  crystal  as  the  two  dimensional  phase  sensitive 
detector.  In  both  drawings  S.M.  denotes  a  scattering  medium. 

Fig  3:  The  experimental  results  for  imaging  an  object  within  a 
scattering  medium:  (a)  The  original  input;  (b)  The  results  after 
transmission  within  the  scattering  medium;  (C)  The  results  using 
a  CCD  camera  as  the  two-dimensional  phase  sensitive  detector. 

Fig  4:  The  same  as  in  Fig  3  respectively  but  for  the  case  when  the 
photorefractive  crystal  has  been  used  as  the  two-dimensional  phase 
sensitive  detector. 
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5.4.2  PATENTS: 


1-  Homodyne  and  heterodyne  imaging  in  a  scattering  medium 

2-  Imaging  in  a  scattering  medium  using  photon  density  waves. 
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